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ABSTRACT 

cn" 

The fluctuations observed in the light curves of some GRB-afterglows (such as GRB 
\ 021004) provide a useful tool to probe the circum-burst density profile and to probe the 

variations in the energy of blast-wave with time. We present a general formalism that 
reduces the calculation of the observed light curve from a Blandford-Mackee blast-wave 
to the evaluation of a one dimensional integral. Using this formalism we obtain a simple 
approximation to the general light curve that arises in more complex situations where 
the afterglow's energy or the external density vary. The solution is valid for spherically 
symmetric profiles and it takes a full consideration of the angular time delay effects. We 
l/-) \ present the light curves of several external density profiles and demonstrate the effects 

of density variations on the light curve. We also re-visit the afterglow of GRB021004 
and we find that the steep decay after the first bump (~ 4000sec) cannot result from a 

o 

of the observed features in the light curve. This may be the first evidence that an angular 
q structure is important in the early stages of the afterglow. 
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spherically symmetric density variation or from the passage of the synchrotron frequency 
through the optical band. This suggests that an angular structure is responsible to some 



Introduction 



The basic theory of the multi wavelength afterglow which follows GRBs is well established. 
The afterglow emission is produced by an interaction between a relativistic expanding fireball and 
the circumburst material. This interaction produces a relativistic blast wave which heats the ex- 
ternal medium and produces the observed emission. The dominant radiation process is most likely 
synchrotron. 

After a short radiative phase the blast wave becomes adiabatic and the hydrodynamic profile 
behind the shock relaxes to a self-similar solution (Blanford & Mackee 1976 ; hereafter BM76). The 
observed synchrotron light curve from a Blanford-Mackee (BM) self similar shell is presented by 
Granot, Piran & Sari (1998) as a two dimensional nontrivial numerical integral. We use here the self 
similar properties of the BM solution to reduce the expression of Granot et. al. (1998) to a simple 
one dimensional integral. The integrand expresses the contribution of the instantaneous emission 
from a BM shell at a given radius 1 . This simplification is important for two reasons. First, it allows 



lr Tlie integrand includes an integration over the angular and the radial dimensions which reduces, due to the 
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an easy calculation of interesting physical quantities (e.g. the afterglow image, see Granot et. al. 
1998) that had to be calculated numerically so far. Second, and more importantly, it enables us to 
obtain a good approximation of the observed afterglow when the external density and/or blast- wave 
energy varies, as long as the variations are spherically symmetric. 

We use our solution to investigate several external density profiles. We find that due to angular 
smoothing, even a short length scale density variation, results in a fluctuation in the light curve 
of more than one order of magnitude in time. Second, we show that above the cooling frequency, 
a density enhancement results in a fluctuation with amplitude of less than 20% , while a density 
drop results in a larger fluctuation (a drop of one order of magnitude in the density results in a 
fluctuation of 40%). 

Nakar, Piran & Granot (2002) investigated different possible explanations to the peculiar after- 
glow of GRB 021004. This unusual afterglow shows a clear deviations from a smooth temporal power 
law decay. A first bump is observed after ~ 4000sec, this bump is followed by a steep decay. Later 
the afterglow shows additional deviations from a power-law decay. Nakar et al. 2002 considered the 
angular effects only approximately. Here we re- visit this afterglow taking a full consideration of the 
angular effects. We show that the first bump and the steep decay that follows it, cannot result from 
spherically symmetric density variations. We show also that the fast decay is inconsistent with a 
simple transition of the typical synchrotron frequency, f m , through the optical band (with a con- 
stant ISM density profile). Although passage of u m is consistent with the timing and with the rising 
phase of the bump (Kobayashi k, Zhang 2002), the fast decay requires another mechanism (most 
likely an angular dependent one). These results suggests that an angular structure is important in 
the early afterglow of GRB 021004. 

In §2 we present our solution. In §3 we consider the effects of density and energy variations. In 
§4 we discuss the afterglow of GRB 021004 and at the last section we present our conclusions. The 
detailed calculations which lead to the solution we present at §2 are described in the appendix. 

2. The Blandford-McKee light curve 

We consider an adiabatic relativistic blast wave propagating into an external medium with a 
regular density profile n ext oc r~ k (k < 3). We calculate the observed flux from a slow cooling 
synchrotron. As commonly done (Sari, Piran & Narayan 1998) we assume that the energy of the 
magnetic field and the energy of the hot electrons, are constant fractions of the internal energy (es 
and e e respectively) and that the hot electrons' initial distribution is a power law with an energy 
index p. We express the solution as a simple one dimensional integral. This solution can be used 
to calculate different properties of the afterglow, or as an approximation to the observed afterglow 
in the cases of a variable external density and/or a variable blast- wave energy (see section 3). 



self-similar structure of the shell, to a general integral that is calculated only once. 
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The spectral shape of the observed flux is well approximated (Sari et al. 1998) by a broken 
power law with several power-law segments. In the slow cooling regime there are four segments 
separated by u a , v m and u c : the self-absorption, synchrotron and cooling frequencies respectively. 
We discuss here only the emission above v a . First we consider the solution when the observed 
frequency, v, is far from the break frequencies and later we address the solution near the break 
frequencies. 

We distinguish between three relative frames. The center of the explosion and the observer 
at infinity are at rest at the source frame. We denote the time in the source frame by t and we 
calibrate t = to the explosion time. The fluid frame is the rest frame of the fluid, we denote all 
variables in the fluid frame with a 'prime'. The observer frame is at rest compared to the source 
frame. However, the observer time is measured according to the arrival time of photons to the 
observer. We denote the observer time as T, and we calibrate T = to the arrival time of a photon 
emitted at t = 0, r = (r is the distance from the explosion). 

We address first the problem of an instantaneous emission from a very thin shell located at a 
radius r and propagating with a Lorentz factor 7 (this problem or related ones where considered 
by Fenimore, Madras & Nayakshin (1996), Kumar & Panaitescu 2000, Ioka & Nakamura 2001, 
Ryde & Petrosian 2002 and others). Due to the relativistic beaming the observer receives mainly 
photons emitted up to an angle of 9 = I/7 relative to the line-of-sight. The observer time delay 
between two photons emitted simultaneously at a radius r, one on the line-of-sight and the other 
from 6 = I/7 is T ang = r/2c^ 2 . Hence, an instantaneous emission in the source frame is observed 
as a finite pulse with a duration ~ T an g . The observed time of the first photon (emitted on 
the line-of-sight) is T[ os = t — r/c, where t is the time in which the shell radiates in the source 
frame. Thus an instantaneous emission from a very thin shell produces a pulse that is observed for 
Tios < T < Ti os + T ang . The emission that arrives at larger T is emitted at larger 9. Fenimore et 
al. (1996) calculated the shape of the pulse for an intrinsic power law spectrum with a power law 
index (3: P v oc v" . They find that the observed flux at T > l) os is: 

F u (T)<x^(l + ±—^) . (1) 

\ 1 ang J 

The observed flux, at an observer time T, from an arbitrary spherically symmetric emitting 
region is given by Granot et. al. (1998): 

FAT) . ^ jf* f rHr £ d(ros ,)=V)f^) J(t - T - — — ), ,2) 

where n' is the emitters density and Pi is the emitted spectral power per emitter, both are measured 
in the fluid frame; 9 is the angle relative to the line of sight, and A -1 = 1/7(1 — v cos9/c) (v is 
the emitting matter bulk velocity) is the blue-shift factor. Below we calculate this expression for 
an instantaneous emission from a shell with a Blandford-Mackee (BM) self-similar profile. The BM 
self similar shell has a width of, AR f» i?/2(4 — k)T 2 , where R and T are the shock's radius and 
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Lorentz factor respectively. BM76 show that all the variables behind the relativistic blast-wave are 
functions only of the dimensionless parameter, X = [1 + 2(4 — k)F 2 ](l — r/ct). x = 1 at the shock 
front and it increases with the distance (down stream) from the shock. 

We will calculate F U {T) by integrating over the contributions from BM shells (Eq. 3 below) 
at different radii. To do so we calculate the contribution to the observed flux at time T of the 
emission emitted by the blast wave during the time that the shock propagates from the radius R 
to R + dR. We denote this contribution as A U (R) • gp(R, T, k)dR where (3 is the spectral index in 
the relevant power-law segment. A u (R)dR is the emitted spectral power during the time that the 
shell propagates from R to R + dR. The function, A v , includes only numerical parameters which 
remain after the integration over r and cosO in Eq. 2 and it depends only on the conditions of the 
shock front along the line-of-sight. The second factor gp is a dimensionless factor that describes the 
observed pulse shape of an instantaneous emission, gp is obtained by integration over cosO and r 
in Eq. 2 and it includes only the radial and angular structure of the shell. The self-similar profile 
of the shell enables us to express gp as a general function that depends only on the dimensionless 
parameter T = T(R, T) (which we define later), hence gp(R, T, k) = gp(T, k). The function gp(T, k) 
is calculated only once for a given (3 and external density profile. 

Using the self similar profile of the emitting region and neglecting terms up to the lowest order 
of 0(7 -2 ) (the detailed derivation is described in details in the Appendix), we reduce Eq. 2 to: 

F ^ = WJ MR)9p(T,k)dR, (3) 

where R ma x satisfies 2 : 

Ti os (R m ax) = T, (4) 

and D is the distance to the source (cosmological factors are not considered through the paper). 
When all the significant emission from the shell at radius R is within the same power-law segment, 
(3, (i.e v is far from the break frequencies) then A v and gp are given by: 



A V (R) = H u { 



R 2 4/3 1/3 ,,-1/3 

ext '<°2+ vl 2 29 sec- cm- Hz' (5) 



where R is the radius of the shock front, n ex t(R) is the external density, E is energy in the blast- 
wave, M(R) the total collected mass up to radius R and H v is a numerical factor which depends 
on the observed power law segment (see Eq. 14). We denote by Q x as the value of the quantity Q 
in units of 10 x (c.g.s). 



2 We denote values at the shock front as functions of the shock radius alone (e.g. Ti OB (R) = Ti os (t(R), R). 
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g(T,p,k) = { 



where 



and 



r l+2(4-fe)T „09,fc) /j 1_ 2(4-fc)f+l \ (2 /3) d < 

(4-fc) Jl X- \ L 2(4-fc) + 2(4-fe) X / X ^ c 

(1 + T)-( 2 -^) 1/ > I/ c 
(T-T los (R)) 



T(R, T) = 



T ang (R) 



li(j3, k) = 3 ■ (71 - 17fc)/(72 - 18fc) - ■ (37 + fc)/(24 - 61fe). 



(6) 



(7) 



(8) 



This set of equations is completed with the following relations between the different variables of the 
blast wave, the observer time and the break frequencies: 



u m = 5- W^nlil.E^M^e^eU Hz, 
v c = 2.5 • W 17 R^nJiy B -2 Hz > 



2 1/2 2 



M(i?) = 47rm p / r 2 n(r)dr 
Jo 

R M(r) 



Ti os (R) 



Jo E{r) V 



3-fc 



-R?~ k 



£(r) 
cRM(R) 



R 



4-fc 



(4- jfe)(3- At)E 
(= 2(4 - fc)T, M (i2)) , 



(9) 
(10) 

(11) 

(12) 

(13) 



2£7(i2) 

where m p is the proton mass and the values in the parenthesis are for a constant energy and a 
density profile of n ex t = n*r~ k (k < 3). The numerical coefficient of Eq. 5, H u , is given by: 



Q-lO-^^eJ/V 1 / 3 v<v n 



2 • 10- 21 (5 • lO 12 )^" 1 )/ 2 e p -_\ e%^ /A ^- p)/2 v m <v<v c ■ 

|±|6 • io~ 3 (2.5 • io 6 f ci 4: 2 2)/4 ^~ p/2 " > ^ 



(14) 



For the two canonical cases of ISM (n ex t = n*, fc = 0) and wind (n ex t = n*r k ; k = 2) Eq. 5 is 



reduced to: 



A u ,i os {R) = H v < 



2 • 10 34 i?i 7 n* o E^ 3 v < 
■W^Rl^nt 3p)/4 EV 2 
(0.07)-p. WRl-^n^'El, v>v c 



(0my?-10™R 2 u 3p n { * 3p)/4 E P 2 v m <v<v c ^—ISM, (15) 

sec ■ cm ■ Hz 



2 • lO 34 ^ 3 n*, 3 4 f < v r 

->\ -p i n 34 p2-p (5-3p)/4 p p 
(0.2)^.10 17 ^n; 2 3 - 3p )/ 4 ^ 2 „>„ c 



V m <V <U C 



erg 



WIND, (16) 
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For v > v c the expression of gp(T,k) in Eq. 6 is analytic. We can obtain an approximated 
analytic expression also for the spectral segment u m < v < v c . The rise time of the pulse is very 
short (g/3 reaches half of its maximal value at T 0.02, and its maximal value at T « 0.14). Hence, 
for T < 0.25 we can approximate gp as a constant. For T > 0.25 the width of the BM shell is 
negligible and we can approximate the pulse decay as an emission from a thin shell with an effective 
angular time, Tang- We calculate Tang by approximating the BM profile as a series of thin shells 
whose Lorentz factors that vary with x, emitting at the same time an radius 3 . Due to the different 
Lorentz factors each shell has its own angular time (T an g(x) °^ X)- The effective angular time is a 
weighted average of T ang (x)- The weights are the emitted spectral power density at \: oc 
(see the appendix for details). Hence: of the 



^(0.25, k) T<0.25 
g/?(0-25,fc)( ^^g ) T>0.25 



T e aU = ^—^T ang {R) v m <v<v c . (18) 

LI — 



The validity of this approximation is shown in Fig 1, which compares the approximation of Eq. 
17 with the complete calculation of Eq. 6 and a full numerical simulation of the emission from a 
BM blast wave. The numerical simulations include both the adiabatic and the radiative cooling of 
the electrons. 



2.1. The light curve in the vicinity of the break frequencies 

The above solution is valid only when the observed frequency is far from any of the break 
frequencies (v m and v c ). To understand the behavior in the vicinity of the break frequencies we 
consider a thin shell with an intrinsic broken power-law spectrum. v' b is the break frequency in the 
shell's rest frame {(3 = j3\ for v' < u' b and (3 = for v' > u' b ). v' h is constant along the shell; at 
the observer frame, however, the break frequency is v = v' b / 'A (A" 1 is the blue-shift factor which 
defined after Eq. 2), which decrease with 6. Hence, it is possible that at first the observed frequency, 
z/, is smaller than v , while at later times it is larger than v . In other words at different times the 
observed frequency is within different power-law segments. In this case Eq. 1 is still valid, only 
P = for u/u h > (1 + (T - T los )/T a ng) and = 02 for vjv h < (1 + (T - T los )/T ang ). 

When the instantaneous emission is from a BM shell then the blue-shift varies within different 
parts of the shell (larger x and/or larger 6 result in a smaller A -1 ). Hence, the observed flux from 
an instantaneous BM shell at a given observer frequency at a given time results from emission in a 
range of fluid frame frequencies. Therefore, it is possible that the flux in a given observed frequency 



3 This approximation is not valid for v < v m . In this spectral range the contribution of shells with large x decay 
more slowly, and the width of the BM shell at late times, but not to late (T « 1) can not be neglected. 
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at a given time corresponds to different power-law segments at the emission from different parts of 
the shell. 

For v <C v c the observed emission arrives from the whole width of the shell. The integration 
over the parameter x within the pulse shape, g, is both along the radial coordinate and the angular 
coordinate (larger x is smaller r and lower 0). This integration depends on the spectral index (3, 
which is different for v < v m = 1/3) and v > u m (/% = (1 —p)/2). The transition from one 
spectral index to the other occur at some critical x value, % m , which satisfies: 

where £ = (37 + fc)/(24 — 6k). Here the pulse shape depends also on Xm and it takes the form: 



A u gp lA (T,k, X m) = j^r, A v<Vm I x -M/fc.*) ( i _ + _ ] dx 




(2-/3; 



v < u c (20) 

The weight of the contribution of each power-law segment is given by the corresponding A v . When- 
ever Xm < 1 than the whole shell emits within the same power-law segment, v > u m , and Eq. 20 is 
reduced to Eq. 6. Similarly, when Xm > 1 + 2(4 — k)T the whole shell emits within the power-law 
segment of v < v m . Eq. 20 provides an exact solution of the spectral break at v m and an exact 
light curve break when v m passes through the observed frequency. Just like Eq. 6, g/3 lt /3 2 (T,k,Xm) 
is calculated only once. However in this case it should be calculated for every T and Xm- 

When v « v c the emission arrives only from a thin part at the front of the BM shell. The 
local spectrum of the emission in the fluid frame vary along the shell, and an exact solution should 
follow the exact profile of the local emission. Hence, there is no simple solution for the exact flux 
in the vicinity of v c (a full solution of the break shape in ISM and wind is presented at Granot & 
Sari 2001). A partial treatment of the break is obtained by taking a sharp transition from v < v c 
to v > v c when v{l + T) = v c , i.e taking the part of v < v c at Eq. 6 for T < v c jv — 1 and the part 
of v > v c for T > v c /v — 1. In this approximation gp 2j p 3 (T, k) is discontinuous (/% = —p/2), but 
the observed light curve and spectral break are rather smooth. 



2.2. The emission from a collimated jet 

So far we have dealt with a spherical symmetric systems. However, in GRBs the relativistic 
outflow is most likely collimated into narrow jets with an opening angle 9j. Our solution is not valid 
if the hydrodynamical parameters depend on the angle from the jet axis. But if they do not, then 
we can easily generalize our results to a jet as long as the observation angle relative to the jet axis is 
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much smaller than 0j. In this case the emission from the edges of the jet at a given R is observed at 
Tj(R,T) = (-f(R)ej) 2 as long as % > 1. Hence, in this case gp (f(R,T) > (*y(R)6j) 2 , fc) = and 
we can use all the above equations with this substitution. The hydrodynamic evolution of such a 
jet is similar to spherical symmetric evolution as long as R satisfies Tj = {^(R)6j) 2 > 1. For larger 
radii the hydrodynamical evolution changes (the jet spreads sideway) and a jet break is observed 
in the light curve. The effect of the cutoff, Tj, on the observed light curve is negligible for v > u m 
and the spherical symmetric solution is valid for any observed time before the jet break time. For 
v < u m the decay is slower, and taking g$ = at the edges of the jet is required also before the 
break. Clearly the whole solution is not valid after the jet break. 



3. Density and Energy variations 

In the previous section we have calculated the observed light curve, for a regular external 
density and a constant energy blast-wave. Consider now the effect of variations in the external 
density or in the energy of the blast-wave. If the variations are not too rapid, then the shell profile 
behind the shock can be approximated by a BM self-similar profile with the instantaneous energy 
and external density. The light curve can be expressed as an integral over the emission from a series 
of instantaneous BM solutions. 

It is worthwhile to explore the conditions in which this approximation is valid. We consider 
first density variations and then we turn to energy variations. When a blast wave at radius R 
propagates into the circumburst medium, the emitting matter behind the shock is replenished 
within AR » i?(2 1 /( 4 " fc ) - 1). This is the length scale over which an external density variation 
relaxes to the BM solution. Our approximation is valid as long as the density variations are on 
a larger length scales than AR . Our approximation is not valid when there is a sharp density 
increase over a range of AR. However, the contribution to the integral from the region on which 
the solution breaks is small (AR/R <C 1) and the overall light curve approximation is acceptable. 
Note, however, that a density jump by more than a factor of ~ 21 can produce a reverse shock (Dai 
& Lu 2002) which breaks the BM profile of the shell and the validity of our approximation. 

A sharp density decrease is more complicated. Here the length scale in which the emitting 
matter behind the shock is replenished could be of the order of R. As an example we consider a 
sharp drop at some radius Rd and a constant density for R > R^- In this case the external density is 
negligible at first, and the hot shell cools by adiabatic expansion. Later the forward shock becomes 
dominant again. Kumar and Panaitescu (2000) show that immediately after the drop the light 
curve is dominated by the emission during the adiabatic cooling. Later the the observed flux is 
dominated by emission from R « Rd, and at the end the new forward shock becomes dominant. 
Our approximation includes the emission before the density drop and the new forward shock after 
the drop, but it ignores the emission during the adiabatic cooling phase. 

A sharp density drop with n ex t oc r~ k with k > 4 or an exponential drop which continues over 
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a long length scale breaks the BM solution and therefore our approximation breaks down. Some 
of these cases can be described by self-similar solutions (Best & Sari 2000 ; Perna & Vietri 2002; 
Wang, Loeb & Waxman 2002). We do not consider these cases here. However our calculations can 
be followed with the new self-similar profiles. 

An additional effect of density variations arises from the relation k = — (R/n ext (R))(dn ext /dR). 
When the density varies so does k. This effect is important only at v < u m . For v > f c , the light 
curve does not depend on k. Fig 2a depicts gp(T,k) when u m < v < u c for different k values. 
It shows that g@ depends weakly on k in this spectral segment. Fig 2b shows that for v < v m a 
re-calculation of k with R is needed. 

Spherically symmetric energy variations are most likely to occur due to refreshed shocks, when 
new inner shells arrive from the source and refresh the blast wave ( Rees & Meszaros 1998, Kumar 
& Piran 2000, Sari Sz Meszaros 2000). Kumar & Piran (2000) show that in such case the solution 
has a smooth transition from the BM solution with the energy of the pre-collision blast-wave ( the 
front shell) to another BM solution with the total energy of the two shells. The collision produces, 
however, a reverse shock whose emission has a lower peak frequency than the forward shock emission. 
Clearly our approximation fails to capture the effect of the reverse shock and it does not capture 
the details of the light curve during the time that the shock crosses the outer shell. 

Our method enables a simple calculation of the observed light curve for a given density and 
energy profile. In Figures 3 and 4 we show the observed light curves for several different density 
profiles (with constant energy). Fig 3 depicts the v m < v < v c light curve for a Gaussian (AR/R = 
0.1) over-dense region in the ISM. Such a density profile may occur in a clumpy environment. The 
emission from a clump is similar to the emission from a spherically over-dense region as long as the 
clump's angular size is much larger than I/7. The different light curves are for a different maximal 
over-densities. We find that a maximal over-density of 5 effects the observed light curve during two 
orders of magnitude in time. The effect of a maximal over-density of 40 is observed during four 
order of magnitude in time (note that in all the cases the width of the Gaussian is similar). Even 
a mild short length-scale, over-dense region (with a maximal over-density of 2) influences the light 
curve for a long duration (mainly due to the angular effects). This duration depends strongly on 
the magnitude of the over-density. Note that due to the nonlinear dependence of the observed flux 
on n ex t a narrower Gaussian (smaller AR/R) with an equivalent amount of mass in the over-dense 
region produces a larger fluctuation. 

Fig 4 shows the observed light curve for several density jumps and drops. In the left panel 
v m < v < v c at all times. In the right panel we compare the effect of similar density jumps on 
the light curve above and below v c . When v m < v < v c there is a transition from one power-law 
to another with the same power law index and a flux ratio factor of (r^/ni) 1 / 2 between the two 
power-laws (as expected according to Sari et al. 1998). The duration of the transition is longer 
for larger density contrasts. This transition is observable for a duration of about two orders of 
magnitude in time for a density contrast of 10 and for one order of magnitude in time for a density 



Fig. 1. — .The general pulse shape in ISM, gp(T,k = 0), for electrons' energy indices p = 2.2 and 
p = 3. The thin solid lines corresponds to Eq. 6, the thick lines are the results of Eq. 17, and the 
squares and the circles are the results of a full numerical simulation of the emission from a BM blast 
wave. Note that the approximation of Eq. 17 is very good during the decay of the pulse. It is almost 
impossible to distinguish between the approximation and the exact calculation of gp{T,k = 0) at 
this phase. At T > 10, the approximation is less accurate, but at this stage the contribution of the 
pulse to the observed light curve is negligible. Similar results are obtained for other k values. 
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Fig. 2. — gp(T,k) for different k values. Left) v m < v < v c \ The differences between the shape of 
the pulses are small. The rise and the decay times are similar for different k values, and the the 
peaks' heights ratios between different k pulses are less then 1.5. Right) v < u m ; The differences 
between the shape of the pulses are large, mainly in the decay time. Therefore, at this spectral 
segment it is important to consider the changes in k when the density varies. 
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contrast of 2. When v c < v (right panel) then after a small fluctuation, the light curve returns to 
behave as if there was no jump (or drop). The width of this fluctuation is between three orders 
of magnitude in time for high density contrast (10) to two orders of magnitude in time for small 
density contrast (2). The maximal amplitude obtained for a density jump is ~ 20% (a deviation of 
~ 0.2mag) while the deviation for a density drop can reach ~ 40% (a deviation of ~ OAmag) 

Now, using Eqs. 3-13 we can approximate the observed light curve for given energy and density 
profiles. However, for a given burst we usually have the observed light curve at hand and not the 
energy and density profiles at the source. In this case we can invert Eqs. 3-13 (numerically) under 
the assumption of either a constant energy or a constant density. Thus, we can find the profile 
of the free variable, which produces the observed light curve 4 . The analytic approximation of gp 
(Eq. 17) greatly simplifies this inversion when v m < v. The observed light curve at a given time 
is a convolution of emission at many different source times (or shock front's radii). Inverting Eqs. 
3-13 requires a de-convolution of the light curve to the emission at different radii. Unfortunately, 
deconvolution amplifies small errors in the observed data and the resulting de-convolved signal (or 
in our case the energy or density profile) is highly sensitive to small variations in the observed light 
curve. 

In some cases the inversion of the observed light curve fails. This usually happens when 
the light curve depicts rapid decay. The angular and radial spreading dictates a fastest possible 
temporal decay (see Eqs. 6 and 17 in which F v oc T~( 2_/3 ) at late times). A faster temporal decay 
is impossible even if the emission from the blast wave completely stops. A faster observed temporal 
decay would result in a failure to invert Eqs. 3-13. This failure implies that a new effect (like 
angular dependence), which we do not consider, must be included. 

4. GRB 021004 

The peculiar afterglow of GRB 021004 was observed on October 4'th 2002. The early optical 
detection (Fox et al. 2002), T ~ 500sec, enabled a detailed observation of this afterglow from a 
very early stage. This unusual afterglow shows clear deviations from a smooth temporal power law 
decay. A first bump is observed at T ~ 4000sec, this bump is followed by a steep decay. Another 
smaller bump is observed at T ~ 7 • 10 4 sec and a possible third one at 3 • 10 5 sec. A steepening 
which may be a jet break is observed at ~ Tdays. Several different mechanisms can lead to these 
observations. Some of the machisms suggested so far are external density variations, angular energy 
structure (patchy shell model), refreshed shocks, and a passage of v m through the optical band 
(Lazzati et al. 2002, Kobayashi & Zhang 2002, Nakar et al. 2002, Holland et al. 2002, Pandey et al. 
2002, Bersier et al. 2002, Schaefer et al. 2002, Heyl & Perna 2002). The last scenario (the passage 



4 This procedure is similar to the one we have used in Nakar et al. (2002), however here we take a complete 
consideration of the angular effects. 
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Fig. 3. — The light curves results from a Gaussian (AR/R = 0.1) over-dense region in the ISM. The 
different thick lines are for a maximal over-densities of 40 (dashed- dot), 5 (dashed) and 2 (solid). 
The thin line is the light curve for a constant ISM density. The inset depicts the ratio of the mass, 
M(R) over the mass of an ISM (without the over-dense region), Mism(R)- 




Fig. 4. — Left) The light curves for several density jumps and drops (n = m, R < Rq; n = r&2, 
R > Ro) for v m <v < u c , Right) The flux deviation from an ISM light curve produced by the same 
density jumps and drops for u m < v < v c (bold lines) and u c < v (thin lines). In both panels the 
time, T, is normalized so Ti os (Rq) = 1. 



13 



of v m ) explains only the first bump, and is combined with the emission of the reverse shock which 
should be dominant till the first bump. 

In the following we apply our method to two possibilities 5 (i): A spherically symmetric (or an 
angular scale larger then I/7) density variations and constant energy; assuming that the optical 
emission is at all time above v m and below v c , (ii) The passage of v m through the optical band, 
assuming an ISM density profile and a constant energy. Figure 5 depicts the best fit that we 
obtained with a spherically symmetric density variations with p = 2. In order to get the fastest 
decay possible after the bump, we stop the emission completely at the peak of the bump. It is clear 
that due to the angular spreading, it is impossible to fit the fast decay after the first bump. The 
fit is even worse with larger p . Figure 6 depicts the best fit of the first bump for a passage of v m 
through the R band in an ISM density profile and with a constant energy Again, due to the angular 
spreading, it is impossible to fit the fast decay We obtain a marginally consistent fast decay only 
if we assume that the emission is completely stopped just when v m is in the R-band. But clearly 
such a coincidence is unlikely. 

The conclusion from these results is that it is unlikely that the light curve of GRB 021004 
results from a spherically symmetric fluctuations. This result provides a new evidence that also in 
the early times of the GRB emission, an angular structure (either of the relativistic wind or the 
circum-burst medium) is involved. 

5. Conclusions 

We have presented a simplified solution of the slow cooling synchrotron emission form a BM 
blast wave. This solution separates the observed flux at a given time to the contributions from 
different BM shells with a different radius of the shock front. Using the self-similar profile of the 
BM shells we have shown that the pulse shape of the emission from BM shells at different radii is 
general (independent of the shock's front radius). We have also presented an analytic expression to 
this pulse shape for v > v m . Thus, this pulse shape could be calculated only once (or the analytic 
expression may be used), and the whole solution turns into a simple one dimensional integral over 
the contributions from different radii. This simplification enables an easy calculation of different 
properties of the afterglow which until now had to be calculate using a complicate and computational 
time consuming simulations. 

The main advantage of our solution is that it enables us to approximate the emission from a 
blast- wave with a varying energy and/or a varying external density, as long as these variations are 
spherically symmetric. The advantage of this solution over the method we have presented in Nakar 
et al. (2002) is the full consideration of the angular effects. 



6 Our spherically symmetric model is not applicable to the patchy shell model and refreshed shocks can not explain 
the fluctuations below the expected power-law decay. 
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Fig. 5. — The best fit to the R-band observations of GRB 021004 (see the references to the GCNs 
at Nakar et al. 2002), using a varying density (and constant energy). We assume that v m < v < v c 
at all times and we take p = 2, E = 10 52 ergs, eb = 0.01 and e e = 0.1. The inset depicts the 
density profile which produces this light curve. Note that at R ~ 1.8 ■ 10 17 The density drops to 
zero (the emission stops). Even though, the resulting light curve cannot fit the steep decay of the 
observations. 




Fig. 6. — The best fit to the first bump in the R-band observations of GRB 021004 (see the references 
to the GCNs at Nakar et al. 2002), using the passage of u m through the optical band in an ISM 
density profile. The different thin lines are for p = 2 (dashed) and p = 2.4 (solid). The thick line is 
for p = 2.4, assuming that the emission is completely stopped just when u m is in the R band. 
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We use our solution to approximate the light curve which results from several density profiles. 
We find out that the duration of fluctuations in the light curve, which results from density variations, 
are long even if the length scale of the density variation, AR, is very short (AR/R <C 1). For 
example a density variation with AR/R = 0.1 and a mild over-density results in a fluctuation 
which is observed for two orders of magnitudes in time. We show also that density variations induce 
mild (~ 30%) fluctuations also above v c . Fluctuation induced by a density drop are larger than the 
fluctuation induced by a density jump. These fluctuations are also observed for about two orders 
of magnitude in time in the case of a sharp density jump, or drop. 

We try to fit the early afterglow of GRB 021004, by a spherically symmetric varying density 
and by the passage of v m through the optical band. Both fits fail to follow the fast decay after the 
first bump in the afterglow. This results suggests that an angular structure within the ejecta or 
within the external density is crucial for the production of the early afterglow of GRB 021004. 



Appendix: Derivation of the Light Curve Formula 

In the appendix we show the details of the calculations which lead to Eqs. 3-14. First we solve 
the problem for v < v c and later we consider the solution for v > v c . 

We start from Eq. 2 which gives the observed flux from an arbitrary spherically symmetric 
emitting region. For cosO = c(t — T)/r we obtain: 

A(r, t) = 7 (1 - v cosQ I c) » (2 7 )^ (l + ^^M) . (21) 

At a given t, the front of the blast wave is at radius R(t) and the emitting region is restricted to 
r < R(t). On the other hand, emission from r < r m i n = c(t — T)= R(t) — c(T — Ti os (R)) would not 
reach the observer at time T. Hence, integrating over cosO [5(cos9) = cS(t — T — rcos9/c)/r], and 
keeping only terms of the lowest order of 7 -2 , Eq. 2 is reduced to: 

o rRmax rR / „ \ P / T — T, (r}\~( 2 ~^ 

"""vl iR L M)K '^ ( 1+ ^r) • (22) 

where R m ax, is the maximal radius of the shock, from which an emission from the blast wave 
contributes to the flux at time T, i.e. Ti os (R max ) = T . We also used the spectrum for v < u c : 
P[,(uA,r,t) = P' v , m (r,t) ■ (i>A/i>' m )P, where v' m is the synchrotron frequency in the fluid frame and 
(3 is the spectral index. 

The dependence of the hydrodynamical parameters in a BM self similar shell on \ 1S (BM76): 
the bulk Lorentz factor, 7(x) = Ty / l/2x = j(R)/y/x, the internal energy density in the fluid rest 
frame, e'(x) = e '(R) " X~^ 17 ~ 4fc )/( 12 -3fc) an d the fluid density behind the shock in the observer frame, 
n (x) = n '(x)l(x) = n(R) ■ ^~( 7 ~ 2fe )/( 4 ~ fc ) ; where 7 (-R), e(R) and n(R) are the hydrodynamical 
parameters values at the shock front (at radius R). Now, We can express the observed synchrotron 
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frequency, u m = 2^v' m (x), and the observed spectral power at this frequency P u ,m = lP' u > m (Sari 
et al. 1998 ) as a function of x- ^m(x) °c 7.67^ oc 7(x) \/ e (x)7m(x) and P v , m {x) °< 7(x) s °c 
7(x)V e (x)) where 7 m is the minimal Lorentz factor of the hot electrons distribution. In the slow 
cooling regime the radiative cooling of the electrons is negligible. The adiabatic cooling of a single 
electron is proportional to e'(x)/ n '(x), hence 7 m (x) = lm{R) ■ x~( 2 + 2/c )/( 12 ~ 3fc ) . Now, we can 
represent all the variables in Eq. 22 as a function of the shock front, R, and the dimensionless 
parameter x (which increase with the distance from the shock front). Integrating over x [dr = 
-dxR/2(4 - k)T 2 ] , using the relation: 

T-T los (r(x)) = Tim X-l m] 
TanMx)) X 2(4-fc) X ' {6) 

{T(R) is defined in Eq. 7) and expressing the density behind the shock as n(R) = 2n ex t{R)T 2 (see 
BM76) we obtain: 

Q C ^"m ax 

F ^ T )~ (4-k)Di J MR)9(s(T(R),k)dR v<v c , (24) 
A V (R) = R 2 n ext (R)P v , m (R) [yj^j) v < v c (25) 



where 



and gp(T(R) : k) is defined by Eq. 6. Note that the emitted spectral power density depends on x 
as: l 2 n'Pl, m {"iv' m )-P oc x -^ k ) is defined in Eq. 8). The values of j{R), is m (R) and P u , m (R) 
could be found for any given external density profile (see Sari et al. 1998 and Nakar et al. 2002) as 
we have done in Eq. 5. 

For v > v c the emitting electrons are cooling fast, and only a very thin layer behind the shock 
contributes to the emission at this spectral regime. Therefore the pulse shape of an instantaneous 
emission from a blast wave at radius R, gp{T, k), is similar to the pulse shape of an instantaneous 
emission from a very thin shell (see Eq. 1) with a spectral index of (3 = —p/2. Finelly, the emitted 
spectral power at the shock front, A V (R), is 

A u (R) = R 2 n ext P^ m uJu^ 2 ^- 1 / 2 v>v c , (26) 

where (3 is the spectral index for v m < v < v c . 
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